Universal features of electron-phonon interactions in atomic wires 
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The effect of electron-phonon interactions in the conductance through metallic atomic wires is 
theoretically analyzed. The proposed model allows to consider an atomic size region electrically and 
mechanically coupled to bulk electrodes. We show that under rather general conditions the features 
due to electron-phonon coupling are described by universal functions of the system transmission 
coefficients. It is predicted that the reduction of the conductance due to electron-phonon coupling 
which is observed close to perfect transmission should evolve into an enhancement at low trans- 
mission. This crossover can be understood in a transparent way as arising from the competition 
between elastic and inelastic processes. 
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INTRODUCTION 

Metallic nanowires, including atomic contacts and 
atomic chains have become ideal systems for testing pre- 
dictions of quantum transport theory Many trans- 
port properties in the normal state like shot-noise 0, 
conductance fluctuations [j| and others in the supercon- 
ducting state p| have been measured with a remarkable 
agreement with the theoretical predictions Q]. A key 
ingredient in this analysis is provided by the set of trans- 
mission coefficients r n (the so-called mesoscopic PIN- 
code), in terms of which a generic transport property 
F(V), where V is the bias voltage, can be written as 
F (V) = E„/(V,t„), /(V,tV being the contribution of 
a single conduction channel The validity of this de- 
composition relies on the usually large difference in en- 
ergy scales between the electronic excitations involved, 
~ eV, and the typical scale W e i for the variation of the 
normal density of states. This one-electron description 
applies as far as charging effects can be neglected due to 
the strong coupling of the electronic states of the atomic 
region with the metallic electrodes. 

Recently there has been much interest in the formation 
of metallic atomic chains using STM and break junction 
techniques |5( allowing the possibility to study electron 
transport through these nearly ideal one-dimensional sys- 
tems. In particular, experiments in Au chains show un- 
ambiguously the presence of features associated to a quasi 
one-dimensional phonon spectrum 6] . The experimental 
conductance exhibits tiny steps at voltages correspond- 
ing to vibrational modes of the finite suspended chain. 
A number of theoretical papers have addressed the ex- 
planation of these basic features 0> 0, 0- It could be 
interesting to investigate whether the type of descrip- 
tion commented above in terms of the system transmis- 
sion coefficients could be valid to analyze the effects of 
electron-phonon (EPH) interactions, which is reasonable 
as the relevant phonon energies Hu ~ 10 meV are still 
small compared to the typical W e i ~ 0.1 — 1 eV in metal- 




FIG. 1: Schematic representation of an atomic wire suspended 
between metallic electrodes. For the theoretical description 
the system is decomposed into L (left), R (right) and C (cen- 
tral) regions. The labels 1 and N denote the outermost sites 
of the wire. 



lie atomic conductors [lOj. Starting from a microscopic 
model and performing reasonable simplifying assump- 
tions we show in this work that the main features due 
to EPH coupling in a single-channel conductor can be 
described by a universal function of its transmission r. 
The sign of the conductance correction, which is negative 
for large transmission, evolves into a positive correction 
for low transmission. This behavior can be understood 
as a competition between elastic and inelastic processes. 
We show that these predictions apply naturally to the 
case of atomic chains of arbitrary length. 



THEORETICAL MODEL 

As a starting point we consider a tight-binding model 
for describing the electronic properties which is well 
adapted to spatially inhomogeneous systems like atomic 
size conductors [lj, H2 • This approach can be straight- 
forwardly extended to incorporate the phonon degrees of 
freedom [jj. We consider the type of geometry schemat- 
ically depicted in Fig. Q consisting of a left and right 
electrode electrically and mechanically coupled to a cen- 
tral atomic size region. The model allows in principle 
to incorporate not only the phonons of the isolated cen- 
tral region (as done in most previous works 0, El) 
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but also to take into account the effect of the electrodes 
in their dynamics. Within this approach the electronic 
Hamiltonian can be written as H e [ — . a hijcl a Cj a , 

where cj CT creates an electron in site i with spin a run 
over sites of the whole system); and equivalently for the 
phonon degrees of freedom we have H p h = J2iPi/2M + 
^2nj UiAijiij, where the atomic displacements Hi and the 
corresponding momenta pi satisfy the quantization re- 
lation [uj,pk] — ihdjk- Both the electronic Hamilto- 
nian and the dynamical matrix A can be decomposed as 
H = Hl + Hr + He + Vlc + Vcr , describing respectively 
the isolated L, R and C regions and the coupling between 
the central region and the electrodes. Finally, the EPH 
coupling H e ^ p h = Tliij \j i^i ~ Uj)Cjdi is obtained by ex- 
panding the electronic Hamiltonian to first order in the 
lattice displacements 0. We shall consider the simplest 
version of this tight-binding model where only the hop- 
ping terms between neighboring sites are non-zero and 
thus the EPH coupling constants |Ay| = dh^jdr are 
non-local and short ranged. 

The transport properties of this model are conveniently 
calculated in terms of non-equilibrium Green function 
techniques. Using the electronic propagators in Keldysh 
space Gff(t,t'), where a, (3 = +,— denote the two 
branches of the Keldysh contour > the current between 
two sites i,j can be written as 



(t,t)-Gjr( t ,t)). 



(1) 



The EPH interaction is included to second order in per- 
turbation theory, which is appropriate in the metallic 
case 0| . The corresponding self-energy insertions can 
be written in terms of the unperturbed propagators as 
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D° k f, with 



denoting the phonon propagators. If the phonons 
are assumed to be in thermal equilibrium [l^j . the fre- 
quency dependent propagators in Keldysh space are re- 
lated to the advanced and retarded ones by D^ a ^ = 

(n B + 5p, + )D°f - (n B + 6 a -)D°j a , where n B is the 
Bose-Einstein distribution. On the other hand, the re- 
tarded and advanced propagators are formally given by 
D°' r ' a (o;) = Uw ± 277) 2 — A] and can be evaluated us- 
ing similar techniques as for their electronic counterparts 
[lfj . Up to this level of approximation, the correction to 
the current due to EPH coupling in the chain geometry 
of Fig. H can be decomposed into elastic and inelastic 
contributions, ic. SI = SI e i + SIi n [15j| . where 
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Tl.r and Jl,r being the tunneling rates coupling the 
chain to the leads (assumed to be energy independent) 
and the corresponding Fermi distributions. Within this 
decomposition the inelastic contribution corresponds to 
the transfer of an electron with the real emission (or ab- 
sorption at finite temperature) of a phonon mode; while 
in the elastic process the interaction with phonons in- 
duces a transition into an intermediate virtual state re- 
sulting in a renormalization of the transmission. 



RESULTS 

Let us first analyze the effect of EPH interactions in 
the conductance for the ideal chain geometry. The main 
features of this case at zero temperature are illustrated 
in Fig. The ratio between the first neighbors EPH 
coupling constant A and the hopping term t within the 
chain is taken to be 7 x 10 -3 while the maximum phonon 
mode of the uncoupled chain is taken as Hujn ~ 0.03i (i.e. 
hiON/t <C 1) in order to reproduce the typical features 
observed experimentally for Au chains 6J. The compar- 
ison with the experimental results is illustrated in Fig. 
13 which will be discussed below. As can be observed 
in Fig. |3 the conductance exhibits small steps which 
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FIG. 2: Features due to EPH coupling in the conductance 
for the ideal chain geometry containing 3 and 4 atoms. The 
insets in the right upper and middle panels show the scaling 
of the total conductance step AG and its width W with the 
number of atoms N. The correlation of these features with 
the phonon spectrum is illustrated by the included average 
phonon density of states within the chain in the lower panels. 
Go denotes the quantum of conductance 2e 2 jh. 



can be associated with the onset of inelastic emission of 
longitudinal phonons. The more pronounced structure is 
associated with the highest phonon mode of the chain but 
smaller features related to lower modes also appear. In 
this zero temperature limit it is possible to correlate the 
width of the steps with the corresponding width of the 
resonances observed in the phonon density of states of the 
chain (see lower panels of Fig. [21 . This width is in turn 
controlled by the mechanical coupling of the chain with 
the bulk electrodes, thus providing and intrinsic width of 
the steps in the conductance 01 ■ Notice that only steps 
associated with phonon modes having a definite parity 
(even or odd) are observed for a given number of atoms 
in the chain. The total step size is found to follow sim- 
ple scaling laws with the EPH coupling A and the highest 
phonon frequency of the chain wjy, i.e SG ~ A 2 /ujn, while 
it increases roughly linearly with N as illustrated in the 
upper rig ht inset of Fig. in agreement with the exper- 
imental observations [6|- Physically, this behavior arises 
from the increase of the interaction probability with in- 
creasing chain length and from the increase of the vibra- 
tion amplitude for decreasing frequency. The width of 
the step is also found to scale simply as ~ 1/N (see inset 
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FIG. 3: (color online) Comparison between the full model 
calculations and the experimental measurements for the con- 
ductance of Au atomic chains. The experimental data corre- 
spond to a typical curve from the experiments described in 
Ref. Q, while the theoretical curve has been obtained for the 
case N = 3 using the set of parameters discussed in the text. 



in right middle panel of Fig. |5J). As can be observed 
in Fig|2l there is a good agreement between our model 
calculations and the typical experimental results of Ref. 
0. The theoretical curve corresponds to a short chain 
(N = 3), and exhibits a small step at low bias in addition 
to the main conductance step. In spite of the experimen- 
tal noise this smaller feature can also be distinguished in 
the experimental curve. 

We have also analyzed deviations from the ideal chain 
geometry, in particular we have studied the effect of a 
zig-zag configuration, taking into account the appearance 
of transversal modes in the chain. We have found that 
these effects do not modify substantially the overall be- 
havior discussed above. The robustness of these results 
and the simplicity of the scaling behavior found for the 
linear chain suggest that the theory could be reduced to a 
minimal model keeping only the essential ingredients. As 
we show below this is still the case when the bare chan- 
nel transmission is reduced with respect to one. Thus, 
assuming that the typical phonon energies are small on 
the scale W e i , the energy dependence of the unperturbed 
electron propagators can be neglected. One can also re- 
late these propagators to the scattering properties of the 
system at the Fermi level by means of the relations 
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(5) 

These coefficients are related to the normal transmission 
t by r = \t LR \ 2 = \t RL \ 2 = 1 - \r L , R \ 2 . Notice that the 
transmission can be reduced either due to asymmetries 
in the coupling to the electrodes or to the presence of 
an impurity atom within the chain. On the other hand, 
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FIG. 4: Schematic representation of a minimal model ac- 
counting for the transition between negative and positive con- 
ductance step as a function of the transmission based in the 
mapping proposed in Refs. [2?^ . A single phonon mode be- 
tween two sites is considered. The two channels correspond 
to the absence or to the presence of an excited phonon. The 
arrows indicate all possible paths for an incident electron: 
elastic and inelastic transmission and reflection. 



if one is interested in the total size of the steps and not 
in their width, the expressions can be greatly simplified 
by taking the limit of negligible broadening of the chain 
phonon states while maintaining their thermal popula- 
tion. In this limit it is possible to factorize the corrections 
to the current into a purely electronic factor and another 
one associated with the coupling to the phonon system. 
Up to a uniform background correction 5Gb, which does 
not exhibit structure at the phonon frequencies, one ob- 
tains from Eqs. llol 111 simple analytical expressions for 
the conductance which for a half-filled band and at zero 
temperature reduce to 

SG el - 5G b = —-T 2 -^Y1 F -p^ 9 d e ^l ~ tor) (6) 



N 



SG m = HLt (1 - t) ^ F ^,r0 (\eV\ - hu; r ) , (7) 



where F P . 



.r)\ 2 . With dj 



denoting the amplitude of the chain atom at site j in 
the mode of frequency uj r (normalization X^i a |r- = 
h/Mw r ). It is interesting to observe that these two 
contributions are universal functions of the transmission 
through the chain. Notice that the inelastic contribution 
is always positive while the elastic one yields always a 
reduction of the total current. On the other hand, while 
the elastic contribution increases monotonously with r, 
the inelastic correction behaves as r(l — r) and thus van- 
ishes in the limit of perfect transmission as in the case of 
shot noise in a quantum point contact |l9| . As a result 
of the competition between the elastic and the inelastic 
parts, the total correction exhibits a change of sign as 
a function of the transmission. In this simplified model 
the step in the current at the phonon frequencies scales 
as t(1 — 2t), i.e. the change in the sign of the correction 
occurs for t = 1/2. This type of behavior can also be 
qualitatively obtained by means of a mapping of the EPH 
problem into a coherent multichannel scattering model, 
as proposed by several authors |2(J . For the simplest case 
of a single localized mode between two atoms in which 
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FIG. 5: (color online) Evolution of the conductance steps as a 
function of the transmission for an atomic chain with N — 3 
and 4. The symbols correspond to the full calculation while 
the full lines were obtained using Eqs. 16171 . 



an incoming electron can be transmitted either elasti- 
cally or inelastically, as schematically depicted in Fig. 0] 
this model would predict 5r e i « — (A/i) 2 r(l — r)(l + 2r), 
where r = \r L . R \ = y/1 - r, and 6r in w (A/t) 2 r(l - r), 
i.e. a change in the sign of the correction for r ~ 0.41; 
while the inelastic reflection is given by Spi n = (X/t) 2 T 2 . 
Thus in the r — > 1 (t — > 0) limit the inelastic tran- 
sitions correspond mainly to back (forward) scattering. 
Although this extremely simple model accounts for the 
overall behavior, the absence of a correct inclusion of the 
restrictions imposed by the Pauli principle is responsible 
for the discrepancies with the correct result (Eqs. (|6I7|) ). 
It is worth noticing that the factor F ep h r in Eqs. l|rJI7|) 
shows explicitly the selection rules of modes with a def- 
inite parity. For the highest frequency mode one has 
a-i,N ~ —CLi+i,N which produces a maximum size for the 
corresponding step. It is also straightforward to show 
that the size of this step increases roughly linearly with 
N. The evolution of the conductance steps as a function 
of the transmission predicted by Eas. l|6l7|l is compared in 
Fig. [5] with the results of the full calculation. As can be 
noticed the behavior of the full calculation is reproduced 
except for the smearing of the steps due to the intrinsic 
width of the phonon modes. 

We have also analyzed in detail the validity of the sim- 
plified model leading to Eas. (|6l7|) for significant devia- 
tions with respect to the half-filled band case. One ob- 
tains corrections to these equations which are of order 
(e/W e i) 2 , where e measures the shift of the band with 
respect to half- filling. We have checked numerically that 
these corrections become of importance only in the limit 
(e/W e i) ~ 1, i.e. close to the nearly empty or full cases. 
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CONCLUSIONS 

We have presented a theoretical analysis of EPH inter- 
actions in atomic size conductors. Our description allows 
to incorporate in an equal footing the effects due to elec- 
trical and mechanical coupling with the electrodes. We 
have shown that in the typical conditions corresponding 
to metallic atomic chains the main features of EPH inter- 
actions can be accounted for by simple analytical expres- 
sions factorized into a purely electronic part and another 
describing the EPH coupling. While the electronic part 
can be written in terms of the bare transmission, the 
phonon part accounts for the selection rules and the de- 
pendence on the length observed for the case of an ideal 
chain. A prediction which deserves further experimental 
investigation concerns the inversion of the conductance 
steps for a single channel at r ~ 1/2. Atomic Au chains 
containing S or O impurities which can reduce substan- 
tially its transmission as well as atomic contacts exhibit- 
ing partially open channels like for Al [l(Tj | are promising 
candidates to test these predictions. 
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